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Abstract. We consider two inverse problems for the multi-channel 
two-dimensional Schrodinger equation at fixed positive energy, i.e. the 
equation — Atp + V(x)ip = Eip at fixed positive E, where V is a matrix- 
valued potential. The first is the Gel'fand inverse problem on a bounded 
domain D at fixed energy and the second is the inverse fixed-energy 
scattering problem on the whole plane R 2 . We present in this paper two 
algorithms which give efficient approximate solutions to these problems: 
in particular, in both cases we show that the potential V is reconstructed 
with Lipschitz stability by these algorithms up to 0{E^^ m ' 2 ^ 2 ) in the 
uniform norm as E — > +00, under the assumptions that V is m-times 
differentiable in L , for m > 3, and has sufficient boundary decay. 



1. Introduction 

We consider the equation 
(1.1) - Aip + V{x)ip = Eip, xeM 2 , E>0, 

where 



1:2 



(1.2) V is a sufficiently regular M n (C)-valued function on 

with sufficient decay at infinity, 

M n (C) is the set of the n x n complex matrices. This equation will also be 
considered on a domain D, where 

(1.3) D is an open bounded domain in R 2 with a C 2 boundary. 

Equation at fixed E can be considered as rather general multi- 

channel Schrodinger (resp. acoustic) equation on D at a fixed energy (resp. 
frequency) related to E. It arises, in particular, as a 2D approximation to 
the following 3D equation 

(1.4) - A xz tl} + v(x,z)tl} = Eip, (x,z) e = D x L, 
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where L = [a,b], a, b £ R, v is a sufficiently regular complex- valued function 
on SI and ip\DxdL = (for example): see |23^ Sec. 2]. In this framework, the 
approximate 2D matrix-valued potential V is given by 

(1.5) Vij(x) = \i5ij + / 4>i(z)v(x,z)4>j(z)dz, x G D, 



for 1 < i,j < n, where n £ N, {^j}j g N is the orthonormal basis of L (L) 
given by the eigenfunctions of — 4^ such that 4>j\dL = 0, ~~^f = ^j&ji f° r 
j € AT, and Sy = 1 if i = j and otherwise. 

In addition, equation (jl.ip can be seen as a particular case of the 2D 
Schrodinger equation in an external Yang-Mills field. 

For equation (jl.ip on D we consider the Dirichlet-to-Neumann map &(E) 
such that 

(1-6) QWMbd) = ^ 

azy 8D 

for all sufficiently regular solution tji of (jl.ip on D = D U where v is the 
outer normal of dD . Here we assume also that 

(1.7) E is not a Dirichlet eigenvalue for the operator — A + V in D. 

This construction gives rise to the following inverse boundary value prob- 
lem on D : 

Problem 1. Given $(E), find V on D. 

On the other hand, for equation (|1.1|) on M 2 , under assumptions (|1.2|) . 
we consider the scattering amplitude / defined as follows: we consider the 
continuous solutions ip + (x, k) of (jl.ip . where A; is a parameter, k £ M 2 , k 2 = 
E, such that 



(1.8) ifj + (x,k) =e ikx I -mV2^e~^f {k,\k\^j 



,i\k\\x\ 



/ 1 \ 

+ o — ;= , as \x\ —> oo, 



for some a priori unknown M n (C)-valued function /, where / is the identity 
matrix. The function / on J£e = {(M) £ I 2 x I 2 : I: 2 = I 2 = J?) arising 
in (jl.8|) is the scattering amplitude for the potential V in the framework of 
equation (jl.ip . 

This construction gives rise to the following inverse scattering problem on 

TO2. 



Problem 2. Given / on find 1/ on I 2 . 
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Problems Q] and [2] can be considered as multi-channel fixed-energy ana- 
logues in dimension d = 2 of inverse problems formulated in |10] in dimen- 
sion d>2. Note that Problems 1 and 2 are not overdetermined, in the sense 
that we consider the reconstruction of a M n (C)-valued function V of two 
variables from M n (C)-valued inverse problem data dependent on two vari- 
ables. In addition, the history of inverse problems for the two-dimensional 
Schrodinger equation at fixed energy goes back to [7] (see also |17| flT] and 
reference therein). Note also that Problem Q] can be considered as a model 
problem for the monochromatic ocean tomography (e.g. see [2j for similar 
problems arising in this tomography). 

As regards efficient algorithms for solving Problems [1] and [2] for the scalar- 
case, i.e. for n = 1, see |16 [ I17 [ [18 } I19j. In addition, as concerns numerical 
implementations of these algorithms for Problem [2] for n = 1, see [H |6], and 
references therein. 

Nevertheless, the fixed-energy global uniqueness for Problem [1] (and for 
Problem [2] with compactly supported V) for n = 1 was completely proved 
only recently in [3]. The reconstruction scheme of [5] is not optimal with 
respect to its stability properties, and, therefore, is not efficient numerically 
in comparison with the aforementioned 2D reconstructions of |16[ [T71 [T8| [19] , 
but it is very efficient for proving some global mathematical results. In 
particular: a related global logarithmic stability estimate for Problem Q] for 
n = 1 was proved in [22j; global uniqueness and reconstruction results for 
Problem Q] for n > 2 were obtained in |23| ; a global logarithmic stability 
estimate for Problem [1] for n > 2 was proved in |24j . In addition, Problem 
[2] with compactly supported V can be reduced, for n > 2, to Problem [TJ as 
in |16| for n = 1. This implies, at least, global uniqueness for Problem [2] 
(in the compactly supported case). On the other hand, the uniqueness for 
Problem [2] fails already for scalar (n = 1) real- valued spherically-symmetric 
potentials V of the Schwartz class on R 2 (see |12|). 

The main purpose of the present work consists in generalizing the afore- 
mentioned reconstruction approach of |18| [19] to the case of Problems [1] and 
[2] for n > 2. As well as for n = 1 this functional analytic approach gives an 
efficient non-linear approximation V appr (x, E) to the unknown V(x) of Prob- 
lems Q] and [2j The reconstruction of V appr (x, E) from $>(E) for Problem Q] 
and from / on for Problem[2]is realized with some Lipschitz stability and 
is based on solving linear integral equations; see Algorithms 1 and 2 of Sec- 
tion 3, Theorems 16. 1| 16.21 and Remarks 16.41 16.51 of Section [6l Among these 
linear integral equations, the most important ones arise from a non-local 
Riemann-Hilbert problem. For the scalar case, Riemann-Hilbert problems 
of such a type go back to [15] , Another important part of these equations 
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is used for transforming &(E) for Problem Q] and / on f° r Problem [2] 
into M n (C)-valued Faddeev function analogues h± on ^e, involved in the 
formulation of the above-mentioned Riemann-Hilbert problem. In addition, 

\\V appr (;E)-V\\ =e{E) 

rapidly decays as E — > +00, where || • || denotes an appropriate norm. In 
particular, e{E) = 0(E~°°) as E — > +00 if || • || is specified as || • ||l°°(£)) and 
V G C°°(R 2 ,M n (C)), suppV C D, for Problem Q] and if || • || is specified as 
II ' IIl°°(ir 2 ) an d V £ 5(M 2 , M n (C)), for Problem [2J where 5 denotes the Sch- 
wartz class. In addition, no reconstruction algorithms for Problems 1 and 2 

- comparable, with respect to their stability, with Algorithms 1 and 2 and 
with an approximation error decaying more rapidly than 0(E~z) as E — > 00 

- are available in the preceding literature, even for V £ C°°(R 2 , M n (C)), 
suppF C C I 2 , when n > 2 (in general). 

In spite of the fact that some excellent properties of Algorithms 1 and 2 are 
proved assuming that V is sufficiently smooth and that E is sufficiently great 
in comparison with (some norm of) V, we expect that these algorithms will 
work rather well even for V with discontinuities and for the case when E is 
not very big in comparison with V. This expectation is based on numerical 
results for Algorithm 2 for the case n = 1; see [B] and references therein. 
Numerical implementations of Algorithm 1 for n > 1 and Algorithm 2 for 
n>2 are in preparation. 

Let us emphasize that in the present work we also develop studies of [23j 
on the 2D multi-channel approach to 3D monochromatic inverse problems for 
equation fjl -4[) . In this connection, the principal advantage of the 2D multi- 
channel Algorithm 1 (see section [3]) in comparison with the 3D algorithm 
of |21j is that Algorithm 1 deals with non-overdetermined data and is only 
based on linear integral equations. High energy error estimates for both cases 
are similar. However, properties of Algorithm 1 of the present work are not 
estimated yet with respect to the approximation level n in the framework of 
3D applications. 

Finally, note that multi-channel inverse problems and their applications 
to inverse problems in greater dimensions were initially considered for the 
one-dimensional multi-channel case, see [I], [26]. As one of the most recent 
result in this direction see |14j . 

Acknowledgements. The first author was partially supported by the 
Russian Federation Government grant No. 2010-220-01-077. We thank V. 
A. Burov, O. D. Rumyantseva, S. N. Sergeev and A. S. Shurup for very useful 
discussions. 
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2. FADDEEV FUNCTIONS 

In this section we recall some preliminary definitions. 

Under assumptions (jl,2p . we consider the Faddeev functions G(x,k) 
e lkx g(x, k), tp(x,k), h(k,l) and related function R(x,y,k) (see [SJ El ESI 
for n = 1): 

(2.2) V(z, fc) = e ifca 7 + / G{x- y, k)V(y)^(y, k)dy, 

i \ 2 r 



(2.3) h{k,l)= — / e" ii:r y(x)V>(2^)^, 

(2.4) R(x,y,k) = G(x-y,k)+ f G(x - Z,k)V(OR(S,y,k)d£ 

where x = {x 1 ,x 2 ),y = G M 2 , A; = (M fo) G C 2 \M 2 , I = (l u l 2 ) G C 2 , 

Imk = Iml ^ and / is the identity matrix. We recall that 

(2.5) (A + k 2 )G(x,k) = 5(x), 

for x G M 2 , k G C 2 \ M 2 , where (5 is the Dirac delta. In addition: formula 
(|2,2p at fixed is considered as an equation for 

(2.6) i>(x,k) = e ikx n(x,k), 

where /i is soug ht in L°°(M 2 ,M n (C)); formula ((231) at fixed fc and 

y is con- 
sidered as an equation for 

(2.7) R(x,y,k)=e ik ^-y\(x,y,k), 

where r is sous ht in Lj 2 oc (R 2 , M n (C)), with the property that \r(x, y,k)\ ->• 
as |x| — )• oo. As a corollary of (12. ip . (12. 2p and (|2.5p . ^ satisfies (jl.ip for 
E = k 2 = k\ + kl and 

(2.8) (A + A; 2 - T/(x))i?(x, y, fc) = <5(x - y), 

for x, y, G M 2 , k £ C \ M 2 . In addition, ft in (|2.3p is a generalised scattering 
amplitude in the complex domain for the potential V. 
For 7 G S 11 = {7 G M 2 : I7I = 1}, we consider 

(2.9) G 7 (x,fc) = G(x,fc + i07), 

(2.10) Ry(x,y,k) = R(x,y,k + iOj), 

(2.11) ijj^(x,k) = e tkx fjLj(x,k), fj, 7 (x,k) = fj,(x,k + iQj), 

(2.12) MM) = M* + *07,/ + i07), 
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where x, y G R 2 , k G R 2 , I G R 2 . 
In addition, the functions 

(2.13) G + (x,k) = G km (x,k) = ~H^\x\\k\), 

(2.14) R+(x,y,k)=R km (x,y,k) 

(2.15) ^+(x,k) = e ikx fi + (x,k), /j,+(x, k) = [i k /\ k \(x, k), 

(2.16) f(k,l) = h km (k,l), 

for x, y,k,l G R 2 , \k\ = \l\, are functions from the classical scattering theory; 
in particular, / is the scattering amplitude of (|1.8p and Hq is the Hankel 
function of the first type. We also define 

(2.17) h ± (k,l)=h ±k± (k,l), 

(2.18) V±(x, k) = n ±k± {x, k), ip±(x,k) = ip ±k± (x,k), 

(2.19) R ± (x,y,k)=R ±k Jx,y,k), 

where k,l,x,y G R 2 , \k\ = \l\, k± = \k\~ 1 (-k 2 , h) for k = (kx,k 2 ). Note 
that /i+ 7^ fi + , ip+ 7^ an d -R+ 7^ -R + in general. We shall consider, in 
particular, the following restriction of the function h: 

(2.20) b(k) = h(k, -k), for k G C 2 , k 2 = E > 0. 
We now introduce the notations 

Z = X\ + iX2, Z = Xl — 

_9__l/_5_ ._d_\ JL-ifJL 9 

dz 2 \dxi 8x2/ ' dz 2 \dxi dx 2 

A = E- l ' 2 {k x + ik 2 ), X' = E-^ih + il 2 ), 

where x = {x u x 2 ) G R 2 , = {h,k 2 ),l = (h,l 2 ) G C 2 , k 2 = I 2 = E G K+. 
In the new notations 

(2.22a) h = i J E 1 /2 (A + A -i ); fc 2 = ^V2( A -i - A), 

(2.22b) /x = -^(A' + A'" 1 ), l 2 = \E^ 2 {^ X - A'), 

(2.22c) exp(ikx) = exp[^£ 1/2 (Az + A -1 z)], 

where A, A' G C \ {0}, z G C and the Schrodinger equation (jl.ip takes the 
form 

8 2 

(2.23) - 4— + V(2)V» = z G C. 
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In addition, the functions / from (|1.8p and (|2.16p . h± from (|2,17p . fj, + ,ip + 
from d2HHJ), fi±,il>± from (l2~18|) . V from (J22J), /i from ([221) and ft from flZZHI) 
take the form 



f = f(X,X',E), 
(2.24) M + = M +(^,A,£;), 



/i ± = /i ± (A,A',£), 

= 4,+ (z,\,E), 



where A, A' G T, z £ C, E £ R+, 

(2.25) n = fj,(z,\,E), if> = ip(z,X,E), b = b(X,E), 
where A G C \ T, zeC,BeK + . Here 

(2.26) T = {C : C e C,|C| = !}• 

Under assumption (jl.2p . for E sufficiently large the function fi(z, A, E) has 
the following properties (see |18| [19] for n = 1 and Section 4 for n > 2): 

(2.27) n(z, A, -E) is continuous in A G C \ T; 

(2.28) //(z, A(l T 0), E) = n±(z, A, E) for AgT; 

(2.29) ^±(z, A, E) = fi + (z, A, E) 

+ 7Ti jf /z + (*, A", £) X+ (±i (^7 " y ) ) ^±(A ; A", z, £)|dA"|, 
for AgT, where 

(2.30) = for s < 0, x+(s) = 1 for s > 0, 



(2.31) h±(X,X',z,E) = exp 



(2.32) 



-^L(j,(z,X,E) = fj, ( 2,-i,£? ) r(A,z,£), 



8X' 



for A G C \ T, where 
(2.33) 

r(A, z, -E) = exp 



A 



£i/2 ( Az + y + Az + 
V A 



A 



sign(AA - l)b(X,E), 



where b is defined by means of (|2.20p and (|2.22ap : 
(2.34) 
(2.35) 



A 



A — > oo, 







V(z) = 2iE 1 / 2 —fi_ 1 (z,E). 
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The following formula is valid (see |19] for n = 1 and Section 4 for n > 2): 

— - y M ( Z ,-^£;)r(C^,E)dReCdImC 

+ 2^ / M^.C^KKI 1 



(2.36) V{z) = HE 1 ! 2 



for z ££, E sufficiently large and L>_ = {C : (eC, |C| > 1}. 



3. Reconstruction algorithms 

We present here Algorithms 1 and 2, which yield approximate but suffi- 
ciently stable solutions to Problems 1 and 2, respectively. These algorithms 
have a final common part: the reconstruction of the approximate potential 
V appr starting from h± of (|2.17p . Thus, for the sake of clarity, we first give 
the different initial parts of the algorithms — that is, the reconstruction of h± 
starting from &(E) for Algorithm 1 and from / for Algorithm 2 — and then 
the final common part. 

Note that in both algorithms we consider in particular the functions 
ip±,h±,fi- of (pUTj) . (p2~18j) and /i + of (pTTSj) . In addition, in Algorithm 
I , in the definitions of these functions we assume that V = on R 2 \ D. 

Algorithm 1 ($>(E) — > h±). Given $>(E), for E sufficiently large, we 
first reconstruct ip±(x, k)\dD, k 6 M. 2 ,k 2 = E, with the help of the following 
Fredholm linear integral equation (see [16] for n = 1 and Section 4 for n > 2): 

(3.1) 4±(x,k)\ dD = e ikx I+ [ A ± (x,y,k)^ ± (y,k)dy, k G R 2 , k 2 = E, 

JdD 

where 

(3.2) A±(x,y,k)= I G?±(a:-e,fc)(*-*o)(e,y,^K, x,y e 3D, 

JdD 

(3.3) G ± (x,k) = G+(x,k)-^- I e iWdx X +(±ek ± )d9, 

4iri Jgi 

I is the identity matrix, (<I> — &q)(x, y, E) is the Schwartz kernel of the oper- 
ator $(E) — <&o(E), &o(E) is the Dirichlet-to-Neumann operator associated 
to the zero potential in D at fixed energy E, G + (x,k) is defined in (|2.13p . 
k± = (—k2,ki) for k = (^1,^2), dy, d£ denote the standard Euclidean mea- 
sure on the boundary dD and dO denotes the standard Euclidean measure 
on S . 
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Then, in order to obtain h±, it is sufficient to use the following formula 
(see |16) for n = 1 and Section 4 for n > 2): 
(3.4) 

h±(k,l) = -—2 / e~ ilx ($-$o)(x,y,E)^±(y,k)dydx, (k,l) € Jt, E - 

{ 27T r JdD JdD 

Algorithm 2 (/ — > h±). Starting from / on (f° r E sufficiently large), 
one directly recovers h± solving the following integral equation (see \17\ [T9] 
for n = 1 and Section 4 for n > 2): 

(3.5) 

h±(X,X',E) - ivi J f(\",A',E) X+ (±i (± - y)) MA, A", E)\d\"\ 

T 

= f(X,X',E), (A, A') G T x T. 

Algorithms 1 and 2 (h± — > V appr ). We begin with the construction of 
jl + , an approximation to fi + of (j2. 15[) ; this is done by solving the following 
integral equation arising from the non-local Riemann-Hilbert problem (|2.27|) - 
(|2.34p for n in the approximation that b = at fixed E (see |19) for n = 1 
and Section 4 for n > 2): 

(3.6) A, E) + / A + (-z, A', E)5(A, A', z, E)|dA'| = I, X e T, zeC, 

where E 1 is sufficiently large and 
(3.7) 

- \ J T h-(cx;,,E )x+ (_,(! - f)) 

-I/ T MC.V.^)x + ( ! (^-f)) fT ^ y , 

where X-h ^± are defined in (|2,30p . (|2,3ip . Then one can obtain an approx- 
imation |Et_ to /i_ via (|2.29p . used as follows: 

(3.8) p.-{z,\,E) = fl+(z,X,E) 

+ fi + (z, A", E) X+ (-i - y) ) MA, A "' ^)l dA "l' 

for A € T, 2 € C. Finally, the approximate potential V appr (-,E) can be 
obtained using the following formula (see [18\ [19] for n = 1 and Section 4 
for n > 2): 

(3.9) V appr (z,E) = 2^ 1/2 ^ (±-,J T ii-(z,(,E)i(\d(\\ . 
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The approximate potential V appr depends in a non- linear way on &(E) 
in Algorithm 1 and on / on ^e in Algorithm 2, in spite of the fact that 
both algorithms are based on solving linear integral equation. In the linear 
approximation near zero potential, the following formulas hold: 

(3.10) 

h ± (k,l)*-±-[ [ e^- lx + k ^(^-^ )(x,y,E)dxdy, (k,l) € Jl E , 
(2iry J 9D Jar, 

for linearised Algorithm 1; 

(3.11) h±(X, A', E) « /(A, A', E), A, A' G T, 
for linearised Algorithm 2; 

(3.12) V appr (z,E) « -E 1 ' 2 [ w(z,X,E)iX\dX\, 

n Jt 

where z £ D for linearised Algorithm 1 and z£C for linearised Algorithm 
2, and 



(3.13) w(z,X,E) = ^LiJ T exp --E 1 ' 2 (Xz + j - X'z - ^ 
x sign (-i " y)) h±(X,X',E)\dX'\j , 

for z G C, A G T, £ > 0. 



3.1. Algorithm 1 with a non-zero background potential A. Consider 
a potential defined as in (|l,5p . where the diagonal matrix A, defined as 
Aij = Aj(5jj, is supposed to be a known background potential. In this case 
Algorithm 1 admits the following effectivisations. 

Let Vi = A on D, V\ = on I 2 \ D. The following parts A and B 
provide two different approaches to the reconstruction of i/j±(x, k)\dD from 
<J>(i£) and of h±(k,l) from ip±(x, fc)|aDi the reconstruction of V appr from h± 
is given after in steps C and D. 

A. &(E) — > h±. Starting from <&(E), for E sufficiently large, we first 
reconstruct ip±(x, k)\o£>, k G M. 2 ,k 2 = E, with the help of the following 
Fredholm linear integral equation (see Section 4): 

(3.14) (Id + (Id - A±)~ 1 5A±)ip±(x, k)\ dD = Vi(x, k)\ gD , 
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where 

(3.15) A 1 ± u(x) = A 1 ± {x,y,k)u(y)dy, x G dD, 

JdD 

(3.16) A 1 ± (x,y,k)= [ G±(x-£,fc)($l-$o)(£,2/,£K, x,y £ 3D, 

JdD 

(3.17) 5A±u(x)= [ G±(x-Z,k)(*i-*)&V,E)u(y)dyde, xedD, 

JdDxdD 

ij)±(x,k)\dD = (Id - A±)~ l (e lkx I) are the functions ip±(x,k)\aD for V = V\, 
($1 — $q)(x, y, E) is the Schwartz kernel of the operator &\(E) — §>q(E), 
($1 — ®){x, y, E) is the Schwartz kernel of the operator &i(E) — $(E), $o(E) 
is the Dirichlet-to-Neumann operator associated to the zero potential in D 
at fixed energy E, &i(E) is the Dirichlet-to-Neumann operator associated 
to the potential V\ in D at fixed energy E and u is a M n (C)-valued test 
function on dD. 

In order to obtain h± we use the following formula (see Section 4): 
(3.18) 

+ TT\2 / e- ite ($i -®Q){x,y,E)8i>±{y,k)dydx, 

for G . ///;. where h±(k,l) is defined as in ([23]), (|2~T7)) with F = Vi, 

8ip±(x, k) = ip±(x, k) — ifi±(x, k) and ip±(x, k) is defined as ip±(x, k) in (|2.2|) . 
flZZHD , (I2J8D with V = Vi. 

B. $(i£) — )• h±. As above, starting from &(E), for I? sufficiently large, we 
first reconstruct ip±(x, k)\dD, k G M 2 , /c 2 = E, with the help of the following 
Fredholm linear integral equation (see [20j for n = 1 and Section 4 for n > 2): 

(3.19) M^^OIao = ^±(x,fe)|9D + / A±(x,y,k)i)±(y,k)dy, 

JdD 

for fc£R 2 , k 2 = E, where 

(3.20) A±(a;,y,A0= / £, fc) ($ - $1) (£, y, x,y G dD, 

JdD 

V>±, i?± are defined as ip±, R± of $TZ§, <^>, ([2~T0l) . (l2~TT]l . (l2~l~8|) . (12/191) 
with V = Vi; (3> ~~ ^i) (a;, y, 2£) is the Schwartz kernel of the operator 
<&(E) — &i(E), &\(E) is the Dirichlet-to-Neumann operator associated to 
the potential V\ in D at fixed energy E. 
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In order to obtain h± we use the following formula (see |20] for n = 1 and 
Section 4 for n > 2): 
(3.21) 

h±(k,l) = h 1 ± (k,l)+—— / 4}Ux,-k,-l)(^-^ 1 )(x,y,E)ip±{y,k)dydx, 



dD J 3D 



for (k,l) G . ///;. where h\.(k,l) is defined as in $T3$), (pTTjl with V = Vi, 
ip^z(x, k,l) is defined as the solution of the following linear integral equation 
(sec |20j for = 1 and Section 4 for n > 2) 

(3.22) i>Ux,k,l) = e ilx I+ [ G T (x- 1 /,fc)y 1 ( I /)Vi(l/,fc,Z)dy, 



where x, k, I G R 2 , k 2 = I 2 > and G T is defined in (|3.3p . 

C. h± — > /2 + . We construct an approximation fl + to /i + of (I2.15P via 
the following integral equation which generalises (|3.6p (see Section 4): 

(3.23) (Id + (Id + B 1 )~ 1 6B)fl + (z, X, E) = ^' + {z, A, E), A G T, z£C, 
for E sufficiently large, where 

(3.24) B 1 u(X) = [ uiX^B^X^^z^^dX'l, 

(3.25) SBu(X) = J u(X')[B(X,X',z,E) - B 1 (X,X',z,E)]\dX , \, 

for A G T, z G C, B 1 (A, A', z, £7) is defined as B(X,X',z,E) in with 
h± = h±, is defined as /x + in f|2 . 6[> . (I2.15P with V = V\ and u is a 

M n (C)-valued test function on T. 

D. fi + — > V appr . The final part of the algorithm is the same as for 
Algorithm 1 with zero background potential. We construct an approximation 
//_ to /i_ using formula (|3.8p and then the approximate potential V appr (z,E) 
via formula (13.91). 



In the linear approximation near the potential V\, the following formulas 
hold: 

(3.26a) 

h±(k,l) ps h l ± (k,l) 

+ 7tW / e- Ux ($-®i)(x,y,EWUy,k)dydx 

( 27r ) J 3D J 3D 

--final I e' ilx ^ 1 -^ )(x,y,E)(ld-A 1 ± )- 1 5A ± ij 1 ± (y,k)dydx, 

(3.26b) 
h±(k,l) w 
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+ 7^2 / / ^ T (x,-k,-l)(<^ - ^i)(x,y,E)i; 1 ± (y,k)dydx, 

{ 27r ) JdD J 3D 

for (k, I) G J(e-, 

(3.27) fi + (z, A, E) « ^' + (z, A, £?) - (Id + fl 1 ) -1 ^ 1 '-^*, A, £), 

for A € T, zGC, 

(3.28) 

A, E)^y}_- (Id + fl 1 )- 1 ^ 1 '" 1 "^, A, E) 
+ niJ T ^+(z, A", E) X+ Li (^7 " y ) ) (h- ~ hl)(X, A", z, E)\d\"\, 
(3.29) 

V appr (z, E)nVi- -E 1 / 2 [ |- ((Id + fl 1 )" 1 * A, £0) »A|dA| 

+ « £l/2 XX^[" 1,+(z - A ''^(-«(^-^)) 

x (/i_ - ^.i)(A, A", z, E)] \d\"\i\\d\\, 
for z £ D and 22 sufficiently large. 

4. Derivation of some formulas and equations of 
Section 2 and 3 for the matrix case 

The following formula and equations will be useful: 

(4.1) ^(x,k) = ip + {x,k) 

+ 2m [ i; + (x,l)5(l 2 - k 2 ) X+ ((l - fc)7)h 7 (Jb,0* J 
Jm. 2 

(4.2) h 7 (k,l) = f(k,l) 

W /<™,0^-*W(™-*h)M^, 



for j £ S 1 , x,k,l £ M?, k 2 = E £ M. + sufficiently large, 

(4.3) %£-(x, k) = -27T / £, ( '-'//(.r. fc + £)H(k, -^(f + 2fc£K, 

(4.4) —=~(k,p) = -27T / ^(fc + z,p + e)H(fc, -e)<5(e 2 + 2^r, 

for j = 1,2, k G C 2 \ M 2 , k 2 = E G M + sufficiently large, x,p G M 2 , where 
H(k,p) = h(k,k — p), S is the Dirac delta and the other functions were 
already defined in Section 2. Formula fj4. 1|) and equations (|4.2|) - (|4.4p are 
proved in El [13] for the scalar case: the proof can be straightforwardly 
generalized to the matrix case, where one only has to pay attention to the 
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order of factors (which is indeed different from the formulation given in the 
quoted papers, but coherent with similar results obtained in |25j). 



Now formula (I2.29P follows directly from (|4.ip . (|2.6p . (|2.31|) using notations 
(|X2"T|) ; equation ([2~32]) follows from (jOJ) taking into account ([2720]) . (pl33|) 
and notations (|2.2ip . In addition, equation ()3.5[) is a direct consequence of 
P~2j) (with 7 = ±/c_l) using notations (pT2"I|) . 

Formula (l2~36l) follows from (l2~35l) . (l2~34l) . (l2~32|) . (f2T28l) . (l2~2Tl) (these can 
be proved exactly as in the scalar case) and the Cauchy-Pompeiu formula 



(4.5) u(A) = — / u(h r - - / t r — , AeD, 



for any sufficiently regular M n (C)-valued it in T>, where dT> is sufficiently 
regular. In addition, formula (|3.9p is just formula (|2.36p without the first 
term in the sum. 

Equation (|3,6p is an approximation of the following (exact) equation for 

(4.6) fi + (z,X,E) + J fi + (z,X',E)B{X,X',z,E)\dX'\=I + ^(z,X,E), 
for A S T, z £ C, where 

(47) (p(z,\,E) — ^(,,_>, S ) r({ ,,,B)^^f. 

The derivation of (|4.6p can be found in |19] for the scalar case and its gener- 
alisation to the matrix case is straightforward (paying attention to the order 
of factors) . 

Formula ()3.3p is a result of [9 J , while formulas (|3.ip , (|3.2p , (|3.4p are results 
of |16| for the scalar case and can be proved for the matrix case following 
the scheme of [23], where similar formulas appear. 



Formulas (I3T4P and (l3~T8|) follows from (EHJ) and (13H) . 

Formulas (|3.19|) - (|3.22|) are results of |20| for the scalar case and can di- 
rectly extended to the matrix case following the scheme of |23| because, in 
particular, the general matrix version of Alessandrini's identity in |23j works 
for our diagonal background potential A. 

Finally, equation (|3,23p follows from 
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5. Function spaces and some estimates 

We introduce some function spaces, which will be useful to prove the 
high-energy convergence of our algorithms. For m G N, e > we consider 

W m ' 1 (M 2 ,M n (C)) = {u : d k u G L 1 (R 2 , M n (C)) for |fc| < m}, 
W™> l {R 2 ,M n {C)) = {u : x £ d k u G L^R 2 , M„(C)) for \k\ < m}, 
(x £ u){x) = (1 + \x\ 2 ) £/2 u{x), fcG(NUO) 2 , \k\ = h + k 2 , 

for a G]0, 1], s£lwe consider 

C a ' s (R 2 ,M n (C)) = {u : ||u|| a> . < oo}, 

where 

ll^ll a,s — || ^ ^||ce 

w(p + ~ w (p)\ 



\H\a= SUp [\w(p)\ + 

P,?GR 2 , lfl<l V I? I 

= (1 + |p| 2 ) s/2 n(p), |u(p)| = max |uy(p)|; 

l<«,j<n 

in addition we consider ^^(R 2 , M n (C)), defined as the closure of C£°(R 2 , M„(C)) 
(the space of infinitely smooth functions with compact support) in || • || as . 
Let 

(5.1) V(p) = / e tpx V(x)dx, p G R 2 . 
If a matrix-valued potential V satisfies 

(5.2) V G W^ X {R 2 , M n (C)) for some e > 0, m G N, 
then 

(5.3) V eH a , s {R 2 ,M n {C)), ae]0,l], s G R+, 
where a = min(l,e), s = m. Let 

(5.4) S(r) = (1 -r)"V. 
We have the following results: 

Proposition 5.1. Lei f/ie condition (|5.3j) 6e ra/id. TTien 

(5.5a) - V"(fc - 0| < E(r)||F|| a , 5 (l + \k- l\ 2 y s / 2 , 

(5.5b) |fl" 7 (fc,p) - V(p)\ < E(r)\\V\\ a Jl +p 2 )- 8 / 2 , 
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forr = \k\- a ci(a,s,a,n)\\V\\ a ,s < 1, k,l,p £l 2 , 7 6 S 1 , k 2 > 1, 
(5.5c) \H(k,p) - V(p)\ < S(r)||F|| a , s (l +p 2 )~ s/2 , 

for r = \Rek\- a Cl (a, s,a,n)\\V\\ a , s < 1, k G C 2 \ R 2 , p G R 2 , R 3 k 2 > 1. 
In particular 

(5.6a) \f(k,l)\ <2\\V\\ atS {l + \k-l\ 2 )- s/2 , kj£R 2 , 

(5.6b) \H y (k,p)\<2\\V\\ atS (l+p 2 )~ s/2 , k,p£R 2 , 7 G S 1 , 

(5.6c) |fl"(fc,p)| < 2||y|| QjS (l+p 2 )- s/2 , k£C 2 \R 2 , p£R 2 , 

for k 2 > Ei = max(l, (2c\ (a, s, a, n)\\V\\ a)S ) 2 ^ a ), where H 1 (k,l) = 
hj(k, k — I), < a < 1, s > 0, < a < min(l, s). 

Lemma 5.2. Under condition (|5.3p . we have the following estimates: 
(5.7a) 

< \k\~ a c 2 (a,s,a,n)\\V\\ atS , 

forx= (x 1 ,x 2 ) G R 2 , k G R 2 , 7 G S 1 , 
(5.7b) 

< |ReA;|~ CT C2(a, s,cr,n)||y|| ajS , 

/or x = (x 1 ,x 2 ) £ R 2 , k £ C 2 \R 2 , k 2 > Ex(a, s, a, n, \\ V\\ a ,s), where 
< q < 1, s > 1, < a < min(l, s - 1). 

Proposition 15.11 and Lemma 15.21 for the scalar case (n = 1) were given in 
|19j and their generalisation to the matrix case (n > 2) is straightforward. 

6. LlPSCHITZ STABILITY AND RAPID CONVERGENCE OF 

Algorithms 1 and 2 for E ->■ +00 

We present here main rigorous results concerning stability and convergence 
of our algorithms in the case of zero background potential for simplicity. In 
addition, we expect that, for potentials of the form (|1.5p . Algorithm 1 with 
non-zero background potential A (see subsection I3.ip will work even better 
than its version with zero background potential. 

Theorem 6.1 (Stability and convergence of Algorithm 1). Let V G W m,1 (R 2 , M, 
m > 3, suppy C D and let $>(E) be the Dirichlet-to-Neumann operator of 
(|1.6p at fixed energy E, where E > E 2 (a, s, o~,n, ||V|| Q)S ), < a < 1, s = m, 
< a < 1 and E is not a Dirichlet eigenvalue of — A + V and — A in D. 
Then V is reconstructed from §{E) with Lipschitz stability via Algorithm 1 
up to 0(E~( m ~ 2 ^ 2 ) in the uniform norm as E — )■ +00. 



\^{x, k) — I\ + 



djlry{X, k) 



+ 



dfi-y(x, k) 



dxi 



\n(x, k) - I\ + 



dfx(x, k) 
dx\ 



+ 



dfi(x, k) 
dx 2 
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Theorem 6.2 (Stability and convergence of Algorithm 2). Let V satisfy 
(|5.2p . form > 3, and let f be the scattering amplitude of (|1.8p at fixed energy 
E > E2(a, s, a, n, || V||q, !S ), where a = min(l, e), s = m and < a < 1. Then 
V is reconstructed from f on with Lipschitz stability via Algorithms 2 
up to 0(E~( m ~ 2 ^ 2 ) in the uniform norm as E — )• +oo. 

The constant E2 of Theorems 16,11 and 16,21 is precisely stated in Remark 
16.31 The Lipschitz stability of Theorems 16.11 and 16 . 2 1 is specified in the proofs 
of these theorems and is summarized in Remarks 16 . 4 1 and 16.51 The error term 
0(E~( m ~ 2 }/ 2 ) of Theorems 16.11 and 16.21 is made explicit in formula (|6.9p . 

Similarly with the presentation of Algorithms 1 and 2 in section 3, we 
separate the proofs of Theorems 16.11 and 16.21 in several steps. 

Proof of Theorem 1 6. 1\ (&(E) — > h±). We have that equation ([3.ip is a Fred- 
holm linear integral equation of second kind for ip±\dD G E 2 (dD), which is 
uniquely solvable with precise data $ — $0 (the proof of the latter fact is the 
same as in the scalar case; see |16|). Therefore the reconstruction of ip± via 
(|3.ip is Lipschitz stable, with respect to small errors in $ — $0 ( m the L 2 
norm of the Schwartz kernel) , 

As a corollary, the reconstruction of h± in L 2 {^e) from <&{E) — $>q(E) 
via equation (|3.ip and formula (|3.4p is also Lipschitz stable. □ 

Proof of Theorem 1 6. 2\ (f — > h±). Estimates (|5.6p and notations (|2.2ip give 

(6.1a) |/(A,A',£)| < 2\\V\\ a , s {l + E\\-\'\ 2 )- s l 2 , A, A' e T, 
(6.1b) ||/|U2 (TxT) < c 3 n\\V\\ a , s E^ 4 , 

for E > Ei, a = min(l,e), s = m. Now, under the assumptions of 
Theorem I6.2[ integral equation (|3.5p is uniquely solvable for h±(\,-,E) £ 
L 2 (T) for A 6 T, E > E\ (this is a consequence of the unique solvability 
of integral equation (|2.2p for E > E\), In addition, by estimate (|6,lbp . 
for E > max(i?i, (-7rc3n||y|| aiS ) 4 ), equation (|3.5p is uniquely solvable for 
h±(\,-,E) G L 2 (T), A G T, and for h±(;-,E) G L 2 (T x T) by the method 
of successive approximations. This implies the Lipschitz stability of the re- 
construction of h± on T x T from / on T x T, with respect to small errors 
in the L 2 norm. □ 

Proof of Theorems 16'. il and 1 6'. 2\ (h± — > V appr ). The proof follows as in the 
scalar case (that was treated in |19j). except for the order of the terms in 
formulas and integral equations. 
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Estimates (|5.6p . formula (|2.16p and notations (|2,2ip give 

(6.2a) \h±(X,X',E)\<2\\V\\ a>s (l + E\X- X'\ 2 )~ s/2 , X, A' € T, 
(6.2b) \\h±\\ L 2(T xT ) < C3n||y|| QjS £ ,_1/4 , 

for E > E±, s = m and 

(6.3) < a < 1 for Theorem 16.11 and a = min(l,e) for Theorem 16.2 
We define the integral operator B(z,E) as 

(6.4) (B{z,E)u)(X) = / u(X')B(X, A', z, E)\d\'\, 



T 



for A € T, where B(X, A' , z, E) is denned in (j3.7|) and u is a test matrix 
function. The following decomposition holds 

(6.5) B(z,E) =C+Q-(z,E) -C-Q+(z,E), 
where 

(6.6) (C^HA^^^L^C, 

(6.7) (Q±«)(A) = Trijf »(A')x+ f±' (^ - j)) ft ± (A,A',2,E)|<iA'|, 

z G C, A £ T, x + , are defined in (|2,30p and (|2.3ip and u is a test matrix 
function. Thanks to (|6,2p . ()6.5p and properties of the Cauchy projectors C± 
(see [19J for more details), B(z,E) satisfies the estimates 

(6.8a) \\B(z,E)u\\ L 2(T) < c 4 n||y|| aiS £;" 1//4 ||'u|| L 2( T) , 

d 



B(z,E)u 
oz 



< C4n\\V\\ a . s E 1/4 ||u||l2 (T ), 

L2(T) 



(6.8b) 

for z £ C, E > Ex, s = m, a as in (|6.3p . 

Now by estimate ()6.8ap . for E > max(£'i, (c4n|| V || as ) ), integral equation 
(|3.6p is uniquely solvable for fl + (z,-,E) G L 2 (T), at fixed z 6 C, by the 
method of successive approximations. This implies the Lipschitz stability of 
the reconstruction of fi + (z, -,E) on T, at fixed z € C. from h± on T x T 
with respect to small errors in the L 2 norm. □ 

Proof of Theorems \6.1\ and 1 6. 2\ (V appr — > V). Our high-energy convergence 
estimate is as follows: 

(6.9) \V(z) - V appr (z, E)\ < c 5 n||F|| a>s £-( s - 2 )/ 2 , 

where z G C, E > E2(a,s,a,n,\\V\\ ayS ), a as in (|6.3p . s = m, 0<<t<1 
(see |19| for complete details). This estimate follows from (|2.29p . (I2.36p . 
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(|3.6p - (|3.9p . (|6.2b|) and the following estimates (whose proofs for n = 1 can 
be found in [19J): 

(6.10) 



oz 


a 




1 

'"C 


E)r((,z 


(6.11) 






E)- 


fi + (z,-,. 


(6.12) 


dz 


■ 


E)- 


- dz (Z ' 



< omWvWcsE-W 2 , 

£ 2 (T,M„(C)) 



for z G C, s = m > 3, E > E?,{<x, s, a, n, \\V\\ a s ), a as in (|6.3h . 



□ 



Remark 6.3. The constant £2 of Theorems 16.11 and 16.21 can be fixed as 
some constant such that E > Ei implies that 



E>E U \n{z,\,E)\<2, 



d . ^ 
-^»{z,\,E) 



< 1, 



\B(z,E)\\% {T) <-, 



d 

jrB(z,E) 
oz 



op 



L2(T) 



1 

5 2' 



for z G C, A G C, where /U and £? are estimated in (|5.7p and (|6.8p . 

Now, let &y t o(x, y, E), x,y G <9-D, denote the Schwartz kernel of the op- 
erator &(E) — &o(E) considered as precise data for Problem 1. Let <E>y 
denote $^0 with some small errors (for the case of Problem 1) and f denote 
/ with some small errors (for the case of Problem 2). Let V^ ppr denote V appr 
reconstructed from <£y via Algorithm 1 (for Problem 1) and from /' via 
Algorithm 2 (for Problem 2). 

The Lipschitz stability of Theorems 16.11 and 16.21 is summarized in the 
following remarks: 

Remark 6.4. Let the assumptions of Theorem 16.11 hold and let 



(6.13) 5 

Then 

(6.14) 



\&y fi (; ; E) - $y )0 (., ■, E) \\ L 2 [dDxdD) < 5 X {V, E, D, n). 



\V(ippr ~ Vappr\\L°°(D) 



< m (V,E,D,n)5. 



Here Si and rji are some positive constants summarizing the Lipschitz sta- 
bility of Algorithm 1. In particular, 



(6.15) 
(6.16) 



5i(y,£,L>,n) > 3%, 
rn.(y,E,D,n)<rRE t 
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as [|$y j o(-, -, E)\\ip fQ£) X Q£)\ — > 0, for some positive (sufficiently small) d® and 
(sufficiently big) rfl, where 5® and rfl are independent of V and E for fixed 
D and n. 

Remark 6.5. Let the assumptions of Theorem 16.21 hold and let 

(6.17) 5=\\f-f\\ LH ^ E) <d 2 (V,E,n). 

Then 

(6-18) e = \\V apP r - Va PP r\\L°°(m?) < V2(V,E,n)5. 

Here 5 2 and 772 are suitable constants summarizing the Lipschitz stability of 
Algorithms 2. In particular, 

(6.19) 5 2 (V,E,n)>5° 2 , 

(6.20) m (V,E,n) < rj$E, 

as \\f\\L 2 (^ E ) ~> 0) f° r some positive (sufficiently small) 5 2 and (sufficiently 
big) r]2, where 5 2 and rfy are independent of V and E for fixed n. 

Note that in Remark 16.51 t ne norm 1 1 • 1 1 L 2 ( j( E ) is identified with || ■ 1 1 L 2 rp x t\ ■ 
The property that H/Hxar^y —> 0, mentioned in Remark 16.51 is fulfilled, 
in particular, for E — > +00, as a consequence of estimate ()6.1bp . On the con- 
trary, the property that ||$y;o(-, •, -E)||l 2 (9DxOD) ~~ ^ 0) mentioned in Remark 
I6.4[ is not fulfilled for E = Ej,j — >■ 00, for any sequence {Ej}j & ?q of positive 
real numbers such that Ej — > +00 as j — > 00, if V ^ 0. In this connection, 
our high-energy conjecture is that 

(6.21) sup||<Iv,o(-,-,£j)llL 2 (0DxaD) < +OO, 

for some {Ej}j & fq dependent on V, where Ej — > +00 as j — > 00. 

Note that the E factor in the right side of (|6.16p and of (|6.2Up is related 
with the choice of the L 2 norm for estimates of the inverse problem data. For 
example, for Algorithm 2, at least in the linear approximation (|3.1ip - (|3.13|) . 
this factor disappear if || • \\l 2 (^ e ) is replaced by || • \\Lf(J? E )i s = m, where 

(6.22) \\u\\ lt ^ e) = sup (1 + £|A-A'| 2 ) S / 2 KA,A')|. 

(A,A')GTxT 
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